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The entropically based nonlinear Adam- 
Gibbs equation is discussed in the context 
of phenomenologies for nonlinear enthalpy 
relaxation within the glass transition tem- 
perature range. In many materials for 
which adequate data are available, the non- 
linear Adam-Gibbs parameters are physi- 
cally reasonable and agree with those 
obtained from linear relaxation data and 
thermodynamic extrapolations. Observed 
correlations between the traditional Tool- 
Narayanaswamy-Moynihan parameters are 
rationalized in terms of the Adam-Gibbs 
primary activation energy (A/i) determin- 



ing how close the kinetic glass transition 
temperature can get to the thermodynamic 
Kauzmann temperature. It is shown that 
increased nonlinearity in the glass transi- 
tion temperature range is associated with 
greater fragility in the liquid/rubber state 
above T,. 
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1. Introduction 



Enthalpy relaxation near the glass transition is nonex- 
ponential and nonlinear. Nonexponential relaxation is 
very common in almost all types of condensed matter, 
but nonlinearity is unusual because for enthalpy relax- 
ation it becomes significant at small departures from 
equilibrium (typically about 2 K in the temperature 
domain). Only nonlinear viscoelasticity in polymers is 
comparably important for such small and practically 
significant perturbations. In this article the phenomenol- 
ogy of nonlinear enthalpy relaxation near T^ is briefly 
reviewed, and nonlinearity is related to other commonly 
observed features of structural relaxation: the nonexpo- 
nentiality just referred to, departures from Arrhenius 
behavior, and the relationship between kinetics and ther- 
modynamics that is the subject of the article by Angell 
in this Special Issue [1]. 

Empirical evidence for nonlinearity is well estab- 
lished for volume relaxation (and thus of enthalpy as 
well since it includes volume). The most direct and 
compelling evidence is the asymmetric approach to 
equilibrium following temperature steps of opposite 



sign, reported by Hara and Suetoshi [2] for a soda-lime 
silicate glass (plotted in Ref. [3]) and by Kovacs for 
poly(vinyl acetate) [4] (presumably independently of 
Hara and Seutoshi given the inaccessibility of Ref. [2]). 
Another indication of nonlinearity was reported earlier 
by Lillie [5], who observed that the isothermal 
Newtonian viscosity of an unstabilized (nonequi- 
librium) inorganic glass changed with (annealing) time 
as the glass relaxed towards equilibrium. By taking the 
time dependence of the viscosity into account, Lillie's 
shear stress relaxation data could be adequately repre- 
sented by the Maxwell relation [3] 



.(O = .(0)exp(^^), 



(1) 



where a = shear stress, 

rjo = Newtonian shear viscosity, 

Goo = limiting high-frequency shear 
modulus. 
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Since G„ for inorganic glasses is a weak function of 
annealing, Eq. (1) implies a time-dependent shear stress 
relaxation time: 



T(f) 



MO 



(2) 



Hopkins [6] observed that for thermorheologically 
simple materials, Eqs. (1) and (2) could be generalized 
to any nonexponential stress relaxation decay function, 
4>it), using the reduced time ^(t): 



at) 



t 

/ 



dt' 

Tin 



(3) 



in which the zero of time is set when the material first 
falls out of equilibrium.' The concept of the reduced 
time has a long history that, although not without inter- 
est, lies outside the scope of the present report. It should 
be mentioned, however, that it has been used by the 
mechanical engineering community for many years, 
who refer to it as an internal time or in terms of stress, 
strain, or material clocks. The need to incorporate and 
quantify nonlinearity in glassy state relaxation phe- 
nomenology was first recognized by Tool [7-9], whose 
work preceded the introduction of nonexponentiality 
and reduced time by several decades. Tool characterized 
the nonequilibrium state of a glass by introducing the 
fictive temperature, 7}, defined as the temperature at 
which the excess nonequilibrium enthalpy (or any other 
property undergoing relaxation) would be the equi- 
librium value. Equilibrium is therefore characterized by 
Tf{t) = T{t). The fictive temperature, and other aspects 
of the phenomenology of nonlinear relaxation, are dis- 
cussed in detail in a recent review article [10] and in a 
book by Scherer [3]. The nonlinear enthalpy relaxation 
analog to Lillie's time-dependent stress relaxation time 
is an enthalpic retardation time, Th, that is a function of 
the time-dependent enthalpy: 



H{t) = Ho + AHcl,{t,TH[Tit), Hit)]] 



(4) 



Such nonlinear behavior can also be linearized using the 
reduced time, a technique first applied to structural 
relaxation by Gardon and Narayanaswamy [11]. For the 



' Starting the reduced time "clock" at a later time can produce artifacts 
arising from the memory effect associated with nonexponential relax- 
ations, in which relaxation from a particular state is determined not 
only by what that state is, but also on how it was reached. For a 
detailed discussion, see Refs. [10] and [13]. 



ubiquitous stretched exponential form of the decay 
function. <p{t,TH) in Eq. (4) is 



<^(0 = exp(-^(0''), 1>;8>0 



where (3 = nonexponentiality parameter. 



(5) 



The nonlinearity of enthalpy relaxation has important 
practical implications. For example, the excess enthalpy 
(and entropy and volume) of a glass shortens the nonlin- 
ear relaxation and retardation times relative to the linear 
(equilibrium) values, resulting in appreciable relaxation 
(physical aging) in the glassy state. If relaxation in the 
glassy state was linear, annealing or physical aging 
would not occur on practically relevant time scales and 
the material engineering of glasses would be much sim- 
pler (and the science less interesting). The enthalpic 
retardation time is also a function of hydrostatic pres- 
sure and, for polymers, of viscoelastic stresses and 
strains (thus the use of stress and strain clocks referred 
to above). As a result, enthalpy relaxation is coupled to 
both the thermodynamic state of the glass and to vis- 
coelastic retardation and relaxation processes. The 
nature of these couplings is complex and currently un- 
der study [12,13], but lies outside the present consider- 
ations and will not be discussed here in any detail. 
However, one important inference can be drawn from 
data on the effect of hydrostatic pressure on enthalpy 
and volume relaxation. Weitz and Wunderlich [14] 
measured the rates of enthalpy and volume relaxation in 
several pressure-densified glasses (PS, aPIVDVIA, K/ 
CaNOs, sucrose and phenolphthalein), formed by cool- 
ing through Tg under various hydrostatic pressures and 
then releasing the pressure. The glass densities at ambi- 
ent pressure increased with cooling pressure (hence the 
name), and the excess enthalpy increased at pressures in 
excess of about 200 IVIPa. The physical aging rate for 
both enthalpy and volume increased with increasing 
density. Similar observations were reported by Prest 
and Roberts [15] for pressure-densified PVC, who 
found that the rate of enthalpy relaxation increased 
with cooling pressure. These results indicate unambigu- 
ously that excess ("free") volume alone cannot be deter- 
mining the retardation time of either volume or enthalpy 
in such glasses. However, an enthalpy- or entropy-de- 
pendent retardation time can readily accommodate the 
data because these quantities include the pressure-de- 
pendent internal energy as well as volume. 

The most frequently used expression for TH[T{t), 
Tf{t)] in Eq. (4) is the Narayanaswamy equation [16] as 
modified by IVIoynihan [17] (henceforth referred to as 
NIVI): 
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(xAh* (\-x)Ah*)\ 



1 > X > 0, (6) 



where A = preexponential constant, 

X = nonlinearity parameter. 

Ah* = effective activation energy just above Tg. 

The combination of Eqs. (3), (5) and (6) with Boltzmann 
superposition is the most frequently used nonlinear phe- 
nomenology for the study of enthalpy relaxation, and is 
referred to here as the Tool-Narayanaswamy-Moynihan 
(TNM) phenomenology. An equivalent but much less 
used expression for Th(0 is the KAHR equation [18] 
(see Ref. [10] for details).^ The TNM parameters for 
about 30 materials have been collected in Ref. [10]. A 
striking feature of these data is the strong correlation 



between the parameters, exhibited in Fig. 1 in the form 
of X vs Ah*, and in Fig. 2 as x vs jS. These correlations 
are somewhat tentative because the TNM parameters 
have large uncertainties compared with those of typical 
linear relaxation parameters, that arise in part because 
of the data analysis, and in part because the uncertain- 
ties themselves are correlated (the fitting parameters are 
not orthogonal in parameter search space). However, it 
can be confidently asserted that the uncertainties are not 
so large as to move a set of parameters from one end of 
the correlation range to the other, without producing fits 
that lie far outside experimental uncertainties in the heat 
capacity data, and the correlations will accordingly be 
accepted here at face value and as reflecting a physical 
reality demanding of explanation. Because of the empir- 
ical and phenomenological character of the TNM 
formalism, however, other expressions based on some 
sort of physical model are needed before any physical 
interpretation can be attempted of the parameter corre- 
lations, or of the parameters themselves. The nonlinear 
form of the entropic Adam-Gibbs (AG) equation is such 
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Fig. 1. Plot of TNM parameters x vs A/i *. The lines are an aid to the eye, and do not 
conform to any specific mathematical expression. 



^ The historical reason for the KAHR formalism not often being used 
is that it is cast in the form of coupled nonlinear partial differential 
equations, whose solution requires more computing time than does 
Boltzmann integration of responses linearized using the reduced time. 
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Fig. 2. Plot of TNM parameters x vs /3. The line is an aid to the eye, and is not a least 
squares fit. 



an alternative, and is the focus of the present paper. In 
view of the finding that the TNM parameters sometimes 
vary with thermal history [10], however, the discussion 
and interpretations proffered here are restricted to relax- 
ation within and very close to the glass transition tem- 
perature range. Parameters obtained solely from physi- 
cal aging data are not discussed. 

2. The Adam-Gibbs Equation 

The original linear form of the Adam-Gibbs (AG) 
equation [19] derives from a generalization of the transi- 
tion state expression 



TflCr) = A exp 






(7) 



where 

Ayu, = transition state activation energy, 

z * = temperature-dependent number of coopera- 
tively rearranging molecular entities, 

^B = Boltzmann's constant. 

The temperature dependence of z * is determined by the 
macroscopic configurational entropy SSX^ [19]: 






A^A 



Sc{T) 



(8) 



where 



entropy of the smallest number of rearrang- 



ing molecular entities, 
A^A = Avogadro's number. 



Equations (7) and (8) yield 



Th (T) = A exp 



B 



TS.iT) 



in which 



B 






(9) 



(10) 



where C = configurational heat capacity at T'2. 
In deriving these expressions the very weak temperature 
dependence of the preexponential factor. A, has been 
neglected'. Equations (7), (8), and (9) express the AG 



The preexponential factor, A , is regarded here to be just as empirical 
as the arbitrary reference temperatures that often replace it in other 
equations in the literature (such as the WLF equation). This position 
is particularly appropriate for the AG equations, since A corresponds 
to the relaxation time in the limit of high temperature, where the AG 
concept of cooperatively rearranging groups can be expected to break 
down. Accordingly, the same symbol for the preexponential factor is 
used in all the equations presented here, despite the fact that its 
numerical value differs for different equations (sometimes substan- 
tially). 
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thesis that, analogous to the Gibbs-DiMarzio theory [20] 
for a thermodynamic glass transition, it is entropy that 
determines the rate of relaxation. The value of S^{T) is 
computed from 



S.{T) 



T 



where 



A Cp(r) 
r 



dr. 



(11) 



ation for the first time [24]. He computed S^{T) from the 
experimentally determined temperature dependence of 
the configurational heat capacity, parameterized in the 
form ACj, = a- bT. Some time later, the present author 
used the hyperbolic form for ACp (see below) to produce 
the simple nonlinear AG equations described next [25]. 
Nonlinear generalizations of the AG equation are 
obtained by making S^ a function of the fictive temper- 
ature, rather than the thermodynamic temperature: 



ACp(r) = configurational heat capacity, 

T2 = temperature at which the configurational 
entropy falls to zero. 

The quantity ACp(r) is usually identified with the ex- 
perimentally observed difference in liquid/rubber and 
glass/crystal heat capacities (but see later). The temper- 
ature T2 is conceptually identical with the thermody- 
namic Kauzmann temperature, T^ (see article by Angell 
[1]), but is distinguished from T^ here because for non- 
linear enthalpy relaxation it is treated as a fitting 
parameter whose numerical equivalence with T^ is to be 
tested experimentally. 

Nonlinear extensions of the AG equation have a his- 
tory that goes back almost to the time of the original AG 
paper. Since these Symposium Proceedings are a cele- 
bration of 40 years of entropy, it is fitting that this 
history be given here, albeit briefly. Plazek and Magill 
[21] observed just one year after AG that, since the 
excess entropy of the glass is independent of tempera- 
ture, the glassy state activation energy could be calcu- 
lated from the parameter B obtained above T^. A ratio of 
glassy state to liquid activation energies of 0.338 was 
predicted, in excellent agreement with the experimental 
value of 0.333. Since this ratio equals the NM parameter 
X [Eq. (6)], it is evident that a link between nonlinearity 
and the AG prediction of Arrhenius behavior below T^ 
was implicit in the very earliest uses of AG. The connec- 
tion was quantified 2 years later by Macedo and 
Napolitano [22], who deduced the temperature depen- 
dence of Sc needed to produce the VTF equation (see 
Ref. [1] and below) and found that the ratio of glassy to 
liquid activation energies was given by (l-^/T^g). This 
result was generalized by Howell et al. [23], who 
derived an expression for the ratio in terms of the fictive 
temperature dependence of S^: 



TuiTJf) =A exp 



B 



TS,{Tt) 



(13) 



The specific form of TH{T,Tf) is determined by the 
temperature dependence of ACp. For many inorganic 
glasses, this is well approximated by the hyperbolic 
expression 



ACp(r) = ^ 



(14) 



The temperature dependence of ACp for polymers is 
weaker, lying between the hyperbolic dependence and 
constancy, but the results discussed here are not partic- 
ularly sensitive to the form of ACp(r) and Eq. (14) will 
be used throughout the present discussion. Insertion of 
Eq. (14) into Eqs. (11) and (13) yields [10, 25-27] 



Th{T, Tf) =A exp 



B 



T{\ - T2/Tf) 



(15) 



Equation (15) is referred to here simply as the nonlinear 
AG equation, despite the objection that other nonlinear 
AG equations result from different functional forms for 
ACp(r) (see below). This nomenclature is countenanced 
because of the convenient fact that the linear from of Eq. 
(15), obtained by placing Tf = T, is the celebrated Vogel- 
Tammann-Fulcher (VTF) equation [28]"* 



To{T) = A exp 



B 



T-To 



(16) 



where the nonlinear enthalpic T2 and linear VTF To 
parameters are physically equivalent but are distin- 
guished here because, as with the relation between T2 
and Tk, it is a matter for experiment to decide whether 
they are numerically equal. An example of an alternative 



ga (glass) 
Zia (liquid) 



= 1 + 



d\nS, 
dTf 



(12) 



Almost 10 years then elapsed before Scherer applied a 
form of the nonlinear AG equation to enthalpy relax- 



* The derivation of the VTF equation from inserting the hyberbolic 
form for ACp(7') into the entropic AG expression (9) has a history that 
also dates back almost to the origin of AG. References are given in 
Ref. [10], but the detailed history is left for another place and time. 
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nonlinear AG equation is that obtained by assuming that 
ACp is independent of temperature: 



To(T, Tf) = A exp 



B 



T\niTf/T2) 



(17) 



Over the narrow temperature ranges near T^ that are 
considered here, Eq. (17) gives fits to experimental data 
that are indistinguishable from Eq. (15), albeit with dif- 
ferent B and T2 parameters. Indeed, the ranges in tem- 
perature are small enough that the nonlinear AG and 
NM equations also give equally good fits. Relationships 
between the NM and AG parameters are readily derived 
by equating the activation energies above and below T^ 
and applying the approximation Tf ~ T ~ T^. 



Ah* = 



din Th(T) 

d(yr) 



and 



._ ( d\nTniT,Tf) \ 



(18) 



(19) 



For Eqs. (6) and (15), this procedure yields [10, 25-27] 
X = - T2/Tf' = - T2/T^ (20) 

and 

Ah* B B 



R 



(1 - TJT,f 



(21) 



where Tf is the glassy state value of Tf, obtained after 
cooling through T^ but before any significant annealing 
or physical aging can occur.^ Equations (20) and (21) are 
found to be accurate when the same enthalpy relaxation 
data are analyzed using the NM and AG equations [10], 
so that AG parameters can be confidently estimated 
from the more widely published NM parameters. 



3. Physical Interpretation of Nonlinearity 
Parameters 

The nonlinear AG equation has physically meaningful 
parameters that allows several questions to be addressed: 

(i) The validity of Eq. (15) can be tested by com- 
paring the best fit nonlinear enthalpic values of T2 with 
available Kauzmann temperatures Tk- 



(ii) Equation (20) reveals that the NM nonlinearity 
parameter, x, is determined by how close the kinetic T^ 
can get to the thermodynamically defined temperature 
of zero excess entropy, T2. This raises the question of 
what determines the ratio T'g/T'2. 

(iii) Equation (10) indicates that Ajx is calculable 
from experimental values of the AG parameter B, if 
ACp(T'2) and Sc* are known or can be estimated. 

(iv) Equations (8), (10) and (14) indicate that z* at 
Tg is given by 



z*(T,) = 



s*Na Rln(W*) Rln(W*) 



5c(r„) Ci\-T2/T,) 



Cx 



(22) 



' This value of Tf is the definition of choice for Tg. 



where W*(>2) = minimum number of configura- 
tions needed for cooperative relaxation. Since large val- 
ues of z*(T'g) imply increased cooperativity, it is of 
interest to compare it with some independent measure 
of cooperativity. 

3.1 Values of T2, To and Tk 

Table 1 summarizes the values of T2 and Tk for the 
handful of materials for which both are reliably known, 
together with VTF values of To obtained from linear 
dielectric and viscoelastic relaxation data above T^. For 
most materials there is extremely good agreement 
between T2 and Tk (well within uncertainties), but there 
is a significant discrepancy for PS. There is also good 
agreement between the nonlinear T2 and linear To 
values, except for PS again. Given the limited number 
and variety of materials it is inappropriate to draw firm 
conclusions from these data, but it is worth noting that 
the best agreement between T2 and Tk is found for B2O3 
and As2Se3, for which the values of Tk are most reliable 
and for which the temperature dependence of ACp is 
very close to the hyperbolic form of Eq. (14). For glyc- 
erol, there is additional excellent agreement between 
the nonlinear AG value of Ti and the value of T2 
obtained from ac calorimetry [30] in which enthalpy 
relaxation is determined in the linear regime of small 
(sinusoidal) temperature perturbations. These observa- 
tions encourage the belief that nonlinearity is indeed 
determined by the ratio of the kinetically determined T^ 
and the thermodynamic Kauzmann temperature Tk- The 
notable exception to this agreement is PS. This dis- 
agreement cannot be ascribed to experimental uncer- 
tainty, because PS has the most reliable value for Tk of 
any polymer, and its enthalpy relaxation parameters 
have been characterized by more independent groups 
than for any other material and found to be in very good 
agreement [10]. A possible reason for Ti being so far 
below Tk for PS is that ACp is not entirely configura- 
tional, as suggested by Goldstein for glasses in general 
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[3 1 , 32] ." If Tk is forced to equal T^ then only about 40 % 
of ACp is configurational, a figure that lies at the low end 
of the range estimated by Goldstein. This interpretation 
is very speculative, however, and should be regarded 
more as a logical and physical possibility than as a 
compelling explanation. 

Table 1. Values of nonlinear AG Ti, linear VTF Tq, and thermody- 
namic T}t 



Material 


T2 


To 


Tk 




(K) 


(K) 


(K) 


PVAc 


225 


238 {s*f 
247 (WLF)' 




aPMMA 


325 


301 (WLF) 
222 (s*) 


335'' 


PS 


210 


323 (WLF) 


270 


BPAPC 


325 


385 (s*) 




B2O3 


321 




335 


AszSej 


237 




236 


Glycerol 


134 
(128 linear)' 


132 (s*) 


135 



" See Ref. [10] for original references to individual values of T2, To 
and Tk- 

^ Obtained from fitting dielectric data to the VTF equation [Eq. (16)]. 
° From the WLF parameters €2^ obtained from fits to viscoelastic shift 
factors, using To = Tg — €2^. For references to the original data, see 
J. D. Ferry, Viscoelastic Properties of Polymers, Third Edition, John 
Wiley and Sons (1980). 

'' The value of Tk for tmcrystallizable atactic PMMA was estimated 
from the value for crystalline syndiotactic PPMA, using the approxi- 
mation that Sc(Tg) is the same for all tacticities and the experimental 
fact that ACp(r) is identical for all tacticities of PMMA [29], so that 
Tg — Tk is independent of tacticity. 
' From ac calorimetry data (Ref. [30]). 



3.2 Values of A/n 

Values of the AG "primary" activation energy, A/x , 
can be computed from Eq. (10) once s^:* and ACp are 
known, although the reliability of such calculations is 
subject to a number of caveats. First, several types of 
molecular entities are involved in a typical cooperative 
relaxation event, and derived values of A/x are some sort 
of average that cannot be expected to have anything 
more than a semiquantitative significance. Second, ACp 
is an extensive variable and the unit of mass must be 
specified. An appropriate unit for polymers is the 



More precisely, relaxation may be determined by only some of the 
degrees of freedom that contribute to ACp. It is also possible that only 
a part of the entropy of fusion is configurational (or determines 
relaxation), in which case values of Tk obtained by assuming that 
100 % of the entropy of fusion is configurational would be too low, 
and the difference between T2 and Tk would be even larger. If the 
configurational fraction of the heat of fusion is the same as the 
configurational fraction of ACp, however, then the original values of 
Tk are clearly valid. 



Wunderlich bead [33], which can be identified with 
segmental units whose relative rotation is a natural can- 
didate for the primary molecular motion. The AG acti- 
vation energy A/x for polymers can then be identified 
with intersegmental rotational energy barriers. The 
question of mass is more problematic for nonpolymeric 
organic and inorganic glasses, and apart from a few 
brief references to literature data the analysis given here 
is restricted to polymers. Finally, a value for Sc * must be 
specified. This has usually been taken to be A:b In 2, 
since there must be a minimum of two configurations 
available for relaxation to occur (those before and after 
rearrangement). For polymers, which are geometrically 
constrained to essentially one dimension, it has been 
argued [10, 25-27] that ^c* is better approximated as 
k^ In 2^. This value is based on the assumptions that a 
minimum of three segments must cooperatively rear- 
range for the classic crankshaft motion to occur, and 
that there are two nonequivalent rotamers per bond. 
This estimate of ^c* is uncertain, but is clearly better 
than k^ In 2 and, given the unknown contribution to 
rotational barriers by interchain interactions and the fact 
that several types of segmental pairs must be averaged 
for a typical polymer, is considered adequate for the 
present discussion. Values of A^jilk^ per bead are sum- 
marized in Table 2 for four polymers, together with the 
number of beads and values of the number of coopera- 
tively rearranging beads at T^{z*, to be discussed later). 



Table 2. Values of A^u Ik^ (i K per bead) for polymers obtained from 
nonlinear B parameter 



Polymer 


Number 


\lJLlh 


, for W* 


z*(Tg) 




of beads' 


In 2 


ln2' 


(beads) 


PVAc 


4 


6.8 


2.3 


40 


PS 


3 


18 


6.0 


13 


aPMMA 


3 


4.9 


1.6 


120 




4 


3.7 


1.2 




BPAPC 


5 


9.8 


3.3 


45 




6 


8.2 


2.7 





' Number of beads per monomer unit. 

The values A^x. calculated for polymers assuming 
^c* = k-R In 2^ are consistent with expected intersegmen- 
tal rotational energy barriers, as modified by interchain 
interactions. The value for PS is rather high, but might 
be associated in some way with constrained motion of 
the bulky phenyl ring as intersegmental rotation occurs. 
The value for A^x. is also physically reasonable for inor- 
ganic silicates, being comparable with the Si - O bond 
energy if ^c* = k-^ In 2 is assumed [24]. For a series of 
lead phosphate and lead/iron phosphate glasses. Sales 
[34] determined the product A^jiSc* and obtained values 
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of Sc* by assuming that Ajx equaled the P-0 bond becomes infinitely broad (;S ^ 0) as z* ^ °°, is not 



energy. For the lead phosphates he obtained a physically 
sensible value for W* = e.xp{s^*/kB) of 4.6 that was 
independent of lead content. This value increased to 24.5 
when 25 % of the lead was replaced with ferric iron, 
consistent with a more cooperative motion being 
required as the geometric constraints caused by the 
(presumably six fold) crystal field stabilized coordina- 
tion of the iron increased. 

3.3 Values of z* at Jg 

Values of the number of polymer beads that coopera- 
tively rearrange at Tg are included in Table 2, and are 
physically reasonable with the exception once again of 
PS. There is a clear correlation between z*(T'g) and the 
stretched exponential parameter jS, which is exhibited in 
Fig. 3 in the form of a plot of jS vs log z *■ The logarith- 
mic scale for z* has no theoretical basis, and is chosen 
simply because it conveniently produces a near linear 
relation. Note that the decay function becomes exponen- 
tial (;S = 1) as z* approaches unity, consistent with 
relaxation being (by definition) no longer cooperative 
when only a single molecular entity is involved. The 
other expected limit, that the relaxation spectrum 



evident although the apparent limit of z* ~ 2000 as 
;S ^ is large enough to be considered macroscopic 
and physically indistinguishable from mathematical in- 
finity. In any event, the clean correlation between z * and 
13 is in accord with physical expectations, and provides 
further evidence that the nonlinear AG parameters are 
physically reasonable. 

3.4 The Ratio T^/T2 and Parameter Correlations 

In seeking the factor(s) that determine how close Tg 
can get to Ti, and therefore the degree of nonlinearity, 
it is natural to consider the primary activation energy, 
Ajx . The hypothesis to be tested is that lower values of 
Ajx allow Tg to more closely approach 7*2. This can be 
assessed by plotting Tg/T2 = (1 - x)"' as a function of 
B ~ x^Ah *. Although the proportionality factor between 
B and Ajx contains the material dependent ratio s^*/ACp 
[Eq. (10)], this ratio cannot be zero and the limit B ^0 
corresponds uniquely to Ajx ^ 0. A plot of T'g/7'2 = 
(1 -x)"' \s B ~ x^Ah* is shown in Fig. 4, from which 
it is clearly apparent that T'g/7'2 does indeed unambigu- 
ously extrapolate to unit as B approaches zero. The 
AG-based hypothesis that Ajx determines rg/7'2, and the 
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Fig. 3. Plot of nonexponentiality parameter fi vs Adam-Gibbs quantity z*. The line is 
an aid to the eye, to suggest that /3— > 1 asz*— > 1, and is not a least squares fit. 
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associated proportionality between B and {TJTx- 1), 
correspond to the inverse correlation between x Ah* 
seen in Fig. 1. This correspondence is readily demon- 
strated once it is recognized that, by binomial 
expansion, the quantity iTJT2 - 1) = (1 -x)"' - 1 is ap- 
proximately equal to (1 - T2lT^ ~ (accurate within the 
typical experimental uncertainty of 10 % in (1 - x)~' for 
0>x> 0.5). There is a suggestion in Fig. 4 that different 
material types fall on different correlation lines, 
although the uncertainties in the parameters make this 
observation less than compelling. Because of these un- 
certainties, and the ambiguity in choosing the appropri- 
ate unit of mass, no attempt is made here to correlate 
either experimental values of ACp or estimates of Sc* 
with the possible different slopes. 

The correlation between the stretched exponential 
parameter (3 and x (Fig. 2) can also be rationalized in 
terms of Tg/T2, if the ancillary assumption is made that 
the large values of z * associated with small values of 
Tg/Ji correspond to an increased degree of cooperativity 
(as suggested by the polymer data plotted in Fig. 3), and 
that low values of (3 are associated with such cooperativ- 
ity. A plot of jS versus Tg/Ji is shown in Fig. 5 for all 
materials for which data are available and, except for the 
network silicate glasses, the correlation is good and 



consistent with the physically reasonable expectation 
(see Sec. 3.3) that jS — > as Tg/T2 — > 1 (a similar trend 
is suggested by some linear dielectric relaxation data 
[35], that are much less uncertain than the nonlinear 
data under discussion here). 



4. Summary 

The entropy based nonlinear Adam-Gibbs account of 
enthalpy relaxation within the glass transition tempera- 
ture range provides a consistent account of the correla- 
tions observed between the TNM relaxation parameters, 
with physically reasonable parameters in most cases. 
Such uniform consistency is perhaps unexpected, since 
the Adam-Gibbs concept of cooperatively relaxing re- 
gions, and the assumption that the size of such regions 
is determined by the macroscopic configurational en- 
tropy, must both be considered more as heuristic aids 
than as rigorous theoretical concepts. The central AG 
results, that nonlinearity is associated with how close 
the kinetic glass temperature Tg is to the thermodynamic 
Kauzmann temperature T^, and that lower values of 
T^/Tk can be identified with increased fragility [36] 
as proposed by Angell [37], suggest that increased 
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nonlinearity should be regarded as a hallmark of greater 
fragility. The correlations between the various nonlinear 
enthalpy relaxation parameters can then be regarded as 
having the same physical cause as those that form the 
basis for the strong/fragile classification scheme for 
linear relaxation phenomena near but above Tg. 
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